洪水影响预报和风险预警理念与业务实践

我国是世界上洪涝灾害频繁而严重的国家之一,洪水预报预警是防汛减灾工作中重要的非工程措施和洪水防御工作的耳目和参谋。从水文行业的视角,回顾了近年来我国洪水预报预警技术与业务进展,分析了当前洪水预报预警工作面临的新形势和新要求,对比国内外同类行业发展查找了存在的差距,阐述了洪水影响预报和风险预警的定义和理念,从顶层对基于影响预报和风险预警的新一代洪水预报预警系统("国家洪水预报预警系统")总体框架进行了研究和设计,一些关键技术成果在大范围洪水早期预警业务实践中得到了探索应用,取得了较好的效果。     

Location and probability of shoal margin collapses in a sandy estuary

Channel bank failure, and collapses of shoal margins and beaches due to flow slides, have been recorded in Dutch estuaries for the past 200 years but have hardly been recognized elsewhere. Current predictions lack forecasting capabilities, because they were validated and calibrated for historic data of cross‐sections in specific systems, allowing local hindcast rather than location and probability forecasting. The objectives of this study were to investigate where on shoal margins collapses typically occur and what shoal margin collapse geometries and volumes are, such that we can predict their occurrence. We identified shoal margin collapses, generally completely submerged, from bathymetry data by analyzing digital elevation models of difference of the Western Scheldt for the period 1959–2015. We used the bathymetry data to determine the conditions for occurrence, specifically to obtain slope height and angle, and applied these variables in a shoal margin collapse predictor. We found 299 collapses along 300 km of shoal margin boundaries over 56 years, meaning that more than five collapses occur on average per year. The average shoal margin collapse body is well approximated by a 1/3 ellipsoid shape, covers on average an area of 34 000 m² and has an average volume of 100 000 m³. Shoal margin collapses occur mainly at locations where shoals take up a proportionally larger area than average in the cross‐section of the entire estuary, and occur most frequently where lateral shoal margin displacement is low. A receiver operating characteristic curve shows that the forecasting method predicts the shoal margin collapse location well. We conclude that the locations of the shoal margin collapses are well predicted by the variation in conditions of the relative slope height and angle within the Western Scheldt, and likely locations are at laterally relatively stable shoal margins. This provides hypotheses aiding the recognition of these features in sandy estuaries worldwide. Copyright © 2018 John Wiley & Sons, Ltd.     

基于FFT-MA的随机介质建模方法

随机介质是描述地球介质小尺度非均匀性的有效模型，在地震散射波场分析、储层描述等领域具有广泛应用，快速准确的随机介质建模方法是开展相关研究的基础与前提.本文首次将FFT-MA（Fast Fourier Transform Moving Average）算法引入到随机介质建模研究中，分析了该方法与传统随机介质建模方法相比具有的优势，并提出了基于FFT-MA的非平稳随机介质建模方法.建模实验表明，与传统的基于谱分解定理的随机介质建模方法相比，基于FFT-MA的方法在空间域产生随机数序列，使随机数序列与结构参数分离，因此，随机数序列与所建模型存在空间上的对应关系，可以分区域建模和局部修改模型.在非平稳随机介质模型建模时，滑动窗口的大小可以根据自相关长度变化而变化，避免了每个采样点都建立一次完整大小的模型，提高了建模效率.因此，FFT-MA随机介质建模方法能准确构建满足自相关函数要求的平稳及非平稳随机介质模型，具有建模效率高、灵活、实用的优点.     

近期长江大通至南京河段潮动力变化趋势与机制

近期长江感潮河段径、潮动力已然发生变化,但其变化机制与趋势有待进一步探讨。通过对长江大通至南京河段的野外调查,并分析了近40年来大通、芜湖和南京站水文资料,探讨了近期该河段的潮动力变化机制与趋势。结果表明:近年来长江大通至南京段潮动力有增强趋势。具体表现为相近径流量条件,潮差平均增大约10 cm,主要分潮振幅增加10%～30%,潮汐形态系数有减小趋势。引起上述变化的原因可能有:(1)近期感潮河段整体冲刷变深导致潮波上溯阻力减小;(2)口外潮汐动力增强以及海平面上升等使潮汐上溯能力增强。此外,长江流域修建了大量水库群,导致该河段径流量变化由自然因素主导变成自然与人为调控共同作用为主,从而影响了潮动力的相对强弱。     

基于Remez迭代算法求解差分系数的双相介质自适应有限差分正演模拟

利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性.     

赤峰中心地震台地面与井下地震观测系统监测能力对比分析

为分析总结地面与井下地震观测系统的特点，以赤峰中心地震台地面与井下观测系统为研究对象，在噪声分析、噪声功率谱分析、地震监测能力和观测动态范围等方面进行分析。结果表明，2套观测系统的RMS均可达到Ⅱ级环境地噪声水平，井下观测系统噪声小于地面观测系统。对2套观测系统的功率谱密度、有效动态范围的对比均表明，井下观测系统的动态范围比地面观测系统超出约10%，因此，井下观测系统地震监测能力优于地面观测系统，井下观测系统能更有效地记录观测数据。     

基于地基GPS-TEC的中国区域地震电离层效应研究

基于中国陆态网络地基GPS-TEC观测，针对2008～2019年发生在中国区域的7个M_S≥6.0地震，采用滑动四分位距法分析了地震前后的电离层扰动时空分布特征。结果显示，5个地震的震前2～6天，GPS-TEC值出现负异常扰动，地震发生期间及震后电离层TEC出现正异常扰动，主要集中在震后2～7天；GPS观测站距离震中越近，垂直上空的TEC扰动越明显，扰动空间最大范围可达2000km。随着震级的增加，震前电离层TEC异常扰动的发生率有所增加，且异常覆盖的范围也有所扩大。因此，认为震前一周内的电离层TEC变化可能提供揭示电离层扰动与地震活动之间关系的线索。     

2020年6月26日新疆于田6.4级地震前中期和短期预测依据与启示

2020年6月26日新疆于田6.4级地震发生在2020年度全国地震危险区内，震前作出了较好的中期（年度）和短期（月尺度和周尺度）预测，是少有的地球物理观测能力较低地区的强震前中期和短临预测较好震例。本文梳理了中期和短期预测的主要依据及其预测效能，研究表明，震前中期异常主要有流动地磁、多方法组合、5.0级地震平静打破、6.0级地震的准周期活动等；短期异常有4.0级地震活动图像、中源地震影响、于田垂直摆倾斜EW、于田GNSS基准站EW位移、和田GNSS基准站EW位移等。在总结震前分析预测过程的基础上，提出针对地球物理观测密度低地区的地震危险区论证和短临跟踪的建议，为该类地区的地震危险区判定及跟踪工作提供宝贵经验。     

Impact of fractional probability distributions on statistics of hydraulic conductivity,dynamics of groundwater flow and solute transport at a low-permeability site

The fractional Brownian motion (fBm) and fractional Lévy motion (fLm) can easily describe the geometry and the statistical structure of hydraulic conductivity (K) for real-world. However, the fBm and fLm models have not been systematically evaluated when building the K field for a low-permeability site. In this study, both the fBm and fLm are used to simulate the low-K field at NingCheGu (NCG), Tianjin, China. Groundwater flow and solute transport are then computed using MODFLOW and MT3DMS, respectively, and the influence of the fBm/fLm models for K on groundwater flow and solute transport is discussed. Results show that the fLm fits better the statistics of the low-K medium than fBm, and the random logarithmic K (LnK) field generated by fLm is more stable because the resultant LnK field captures more of the measured properties at the field site than that generated by fBm. In contrast, the LnK generated by fBm is more likely to form both high-K channels and low-K barriers. The fBm therefore predicts more extreme behaviours in flow and transport, including the preferential flow, low-concentration blocks and solute retention. The overall groundwater renewal period and solute travel time for the fLm simulation are slightly shorter than those for fBm. The impacts of the fLm and fBm models on the statistics of the resultant LnK fields and the dynamics of groundwater flow and solute transport revealed by this study shed light on the selection and evaluation of the fractional probability distribution models in capturing the K fields for low-K media.     

VARIATION OF THE TRAVEL TIME DIFFERENCE BETWEEN TELESEISMIC PS CONVERTED WAVE AND PP TRANSMITTED WAVE IN FOCAL REGION BEFORE AND AFTER WENCHUAN MS8.0 EARTHQUAKE

When P waves from distant earthquakes meet a velocity discontinuity in the earth's crust and upper mantle, they give rise to a series of converted PS waves besides PP refracted waves. It is possible to monitor the variation of the physical properties of the medium in the limited formation space above the transition zone of the seismogenic zone by measuring the time difference between the teleseismic PS converted wave and the first arrival P wave, that is, time-variation Δt_PS=t_PS-t_P. The advantage of this method is that the transition point of the teleseismic source with similar source is relatively stable at the transition interface, and the accuracy of the measured relative time is high, and the change of the medium in a small range of the seismogenic zone above the conversion interface can be monitored. This paper studies the variation of the travel time difference Δt_PS in focal region before and after Wenchuan M_S8.0 earthquake. We select 2001 to 2012 as the research period, use teleseismic waveforms which occurred in the southern region of Sumatra and Hindu Kush area recorded by Sichuan station YZP and JJS. These teleseisms satisfy 5.0 ≤ M ≤ 6.5, and their waveform signal-to-noise ratio is high with clear initial P-wave motion. The epicentral distance of teleseisms is less than 3 degrees. Then we obtain the variation of the travel time difference Δt_PS between teleseismic PS converted wave and PP transmitted wave recorded during the study period in the two stations. The results show that there is a slow increasing trend of Δt_PS before 2006, and an obvious low value process of Δt_PS appeared in the period about 2 years before the Wenchuan earthquake. The maximum decline was about 0.2~0.3s, more than 4~5 times the measurement error. The low value has a certain degree of return about 2~3 months before the earthquake. The change of arrival time difference indicates that the medium is in different states in different periods of seismogenic process. The sharp decrease in Δt_PS from 2006 to January 2008 may be due to the strong disturbance caused by the stress accumulation of the medium. At this stage, the velocity of P wave and S wave increases with the increase of stress, and the increase of S wave velocity will result in the decrease of Δt_PS. The change of Δt_PS is greatly affected by S wave velocity, so Δt_PS appears to decrease rapidly. Regarding the low value that has a certain degree of return about 2~3 months before the Wenchuan earthquake, the possible reason is that the release of stress is much higher than the accumulation of stress in meta-instability stage. At this stage, the velocity of S wave decreases and the decrease of S wave causes Δt_PS to increase. Then, the Wenchuan earthquake of magnitude 8.0 occurred. It is shown that the teleseismic converted wave method in this paper can monitor the variation of medium's wave velocity before large earthquakes, and it has a good prospect in seismic monitoring and worth further experimental study.